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Recently Jeong and Kim [Phys. Rev. E 66, 051605 (2002)] investigated the scaling properties of 
equilibrium self-flattening surfaces subject to a restricted curvature constraint. In one dimension 
(ID), they found numerically that the stationary roughness exponent a « 0.561 and the window 
exponent S ~ 0.423. We present an analytic argument for general self-flattening surfaces in D 
dimensions, leading to a = Doq/(D + oq) and 5 — D/(D + oiq) where cto is the roughness exponent 
for equilibrium surfaces without the self-flattening mechanism. In case of surfaces subject to a 
restricted curvature constraint, it is known exactly that Qo = 3/2 in ID, which leads to a — 3/5 
and 8 = 2/5. Small discrepancies between our analytic values and their numerical values may be 
attributed to finite size effects. 
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Fluctuation properties of equilibrium surfaces have 
been studied extensively for the last few decades [1] . Sur- 
face roughness is well documented and classified into a 
few universality classes. The Edwards- Wilkinson (EW) 
class is generic and robust for equilibrium surfaces with 
local surface tension [2]. The EW surfaces can be de- 
scribed by the continuum Langevin-type equation 

MM = - v V*h(F,t) + r,(f,t), (1) 

where h(f, t) is the height at site r and time t, rj(f, i) 
is an uncorrelated Gaussian noise, and v represents the 
strength of local surface tension. 

The surface fluctuation width W(L,t), defined as the 
standard deviation of the surface height h(f, t) starting 
from a flat surface of lateral size L, satisfies the dynamic 
scaling relation 

W(L,t) = L a f(t/L z ), (2) 

where the scaling function f(x) — ^ const, for x 3> 1 and 
f(x) ~ x 13 (J3 = a/z) for x < 1 [1,3]. 

The EW universality class is characterized by the dy- 
namic exponent z = 2 and the roughness exponent 
a = (2—D)/2 (D: substrate dimension). In the absence 
of local surface tension (y = 0) , higher-order local sup- 
pression terms like V 2m h (m = 2, 3, . . .) become relevant 
to determine the scaling properties of surface roughness. 
In this case, the scaling exponents become z = 2m and 
a = (2m - D)/2. 

Recently Kim, Yoon, and Park [4] introduced a global- 
type suppression (self-flattening) mechanism which re- 
duces growing (eroding) probability only at the globally 
highest (lowest) point on the surface. They found that 
this global-type suppression changes the scaling proper- 
ties of the EW surfaces: z « 3/2 and a = 1/3 in ID and 
z w 5/2 and a = 0(log) in 2D. The ID roughness expo- 
nent a = 1/3 characterizing the stationary surface fluctu- 
ations could be derived exactly by mapping the surface 
evolution model onto the static self-attracting random 
walk model [5,6]. 



More recently Jeong and Kim (JK) [7] investigated the 
effect of the self-flattening mechanism on the so-called 
restricted curvature (RC) model [8]. The RC model is 
known to have no local surface tension term (y = 0) 
and its dominant suppression term is of the fourth or- 
der (V 4 /i). Accordingly, the ordinary RC model have the 
scaling exponents z = 4 and a = (4 — D)/2. Using Monte 
Carlo simulations, JK found for the ID self-flattening RC 
model that 

z = 1.69(5), a = 0.561(5), and (3 = 0.332(5). (3) 

Again, the self-flattening mechanism changes the scaling 
properties of the RC surfaces. 

JK also studied the height-height correlation function 
G(r), defined as the average of the square of height dif- 
ferences at two sites separated by a distance r. They 
found an extra length scale £ (window size) where the 
correlation function starts to saturate. The correlation 
function in the steady state satisfies the scaling relation 

G(r) = L 2a g (r/£) with £ - L 5 , (4) 

where the scaling function g(x) — » const, for x 3> 1 and 
f(x) - x 2a '-^ (a' = a/6) for x < 1. This type of a 
crossover scaling, due to the existence of a smaller length 
scale (6 < 1) compared to the system size, has been pre- 
viously identified in the so-called even-visiting random 
walk (EVRW) models (see Eq. (29) in reference [9]). In 
fact, the EVRW model is intimately related to the EW- 
type surface model with the self-flattening mechanism [4] . 
JK's numerical estimates for the exponents are 

5 w 0.423, 0.868, and a' = a/6 w 1.33. (5) 

In this Comment, we present an analytic argument 
that predicts the values of the scaling exponents, a and 
6, associated with the stationary properties of surface 
fluctuations. We consider a general equilibrium surface 
growth model, of which the stationary roughness expo- 
nent is known exactly as ao- The partition function for 
equilibrium self-flattening surfaces of lateral size L can 
be written as [4] 
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where the summation is over all possible height configu- 
rations C subject to a given constraint, K is a temper- 
aturelike parameter, and H(C) is the height excursion 
width (the globally maximum height minus the globally 
minimum height) for a given configuration C. Global sup- 
pression for self-flattening dynamics is simply Metropolis 
type evolution algorithm with this partition function to 
reach the equilibrium [4]. For the EW surfaces, one can 
take the surface height configurations subject to the re- 
stricted solid-on-solid (RSOS) constraint, where the step 
heights are allowed to take finite values. In the case of 
the RC surfaces, the local curvature V 2 h is restricted to 
be finite. 

One can decompose the configurational space into sec- 
tors with a constant excursion width H . Then, the par- 
tition function can be rewritten as 
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where u>l(H)cIH is the number of configurations with the 
height excursion width between H and H + dH. One can 
define Ql(H) as the number of configurations with the 
height excursion width less than H as 
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dH' uj(H'), 



(8) 



and it is clear that fi^oo) = Zl(K = 0). 

We estimate SIl(H) as follows. Consider a flat surface 
between two walls separated by a distance H and parallel 
to the surface. The surface starts to fluctuate, following 
its ordinary evolution dynamics. Whenever the surface 
hits either of the two walls, the number of possible con- 
figurations inside the walls will be reduced by a constant 
factor, compared to the no- wall case (H = oo). This 
entropic reduction can be roughly translated as 



n L (H) w fii(oo) cxp[-aA c 



(9) 



where N c is the typical number of contacts between the 
walls and the surface in the stationary state and a is a 
positive constant of 0(1). Typically, there will be 0(1) 
contacts over a block of lateral size £, within which the 
stationary surface width Wo (£) ~ £ a ° is of the same order 
of magnitude as the wall distance H, i.e., 



N c ~ (L/£) D ~ L D /H D / aa . 



(10) 



With the above estimates, the partition function be- 
comes for nonzero K 

Z L {K) = Z L {Q)j™ dH (J 5 e- aLD l HDIao ^ e~ KH , 



= Z L (0)K f 
Jo 



dH e -KH-aL°,H°>~ r 



(11) 



This integral can be evaluated by the saddle point 
method in the limit of large L. The maximum contri- 
bution comes from 



H* 



j^D \ aa/{D+a a ) 



(12) 



The stationary surface fluctuation width should be pro- 
portional to this typical height excursion width, so we 
predict that 



W — K- a o/(D+a ) L a 



(13) 



with a — Dao/(D + ao). 

The length scale I arising naturally in our argument 
should be proportional to the window size £ in the cor- 
relation function. Inside of this length scale, the surface 
does not feel the global self-flattening suppression. As 
our length scale I is explicitly related to the wall distance 
in Eq. (10), we can also predict that 



(14) 



with S = D/(D + a ). 

In summary, for the general self-flattening equilibrium 
surfaces, we predict that 



a 



Dap 
D + ao 



D 



a 



, and a' = - = a , (15) 
V + ao o 



where a is the stationary roughness exponent for the 
equilibrium surfaces without the self-flattening mecha- 
nism. 

For the EW surfaces, a = 1/2 in ID, which leads 
to a = 1/3 and 6 = 2/3. These results agree with those 
by the other analytic (healing time) arguments [9,10] and 
those by numerical simulations [11]. For the RC surfaces, 
ao = 3/2 in ID, which leads to a = 3/5 and S = 2/5. 
The short range behavior governing the value of k should 
be identical to the ordinary RC model, so we also ex- 
pect that k — 1. These results have small discrepancies 
from the JK's numerical results in Eqs. (3) and (5). We 
believe that this may be due to rather small sizes used 
in their numerical simulations. The RC surfaces have 
a Q = (4-D)/2 in D dimensions, so a = D(4-D)/(4+D) 
and S = 2D/(4 + D). It may be interesting to check these 
predictions for the RC self-flattening surfaces in D = 2 
and 3. 

At the upper critical dimensions (D = 2 for the EW 
and D — 4 for the RC surfaces), the surface roughness be- 
comes logarithmic (ao = 0) and the self-flattening mech- 
anism induces only corrections to scaling in the surface 
fluctuations. The dominant correction seems to be inde- 
pendent of system size L [4], but its functional depen- 
dence on K is not fully explored. In higher dimensions, 
the EW (D > 2)and the RC (D > 4) surfaces are asymp- 
totically flat (ao < 0). The self-flattening mechanism 
induces power-law-type corrections to scaling as given in 
Eq. (15). 
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We could not present any analytic explanation for the 
dynamic exponents z and (3. This may be due to the lack 
of a continuum-type equation to govern the self- flattening 
mechanism. It would be very interesting to find such an 
equation which should contain a global coupling term in 
space. 

We thank J. M. Kim and H.-C. Jeong for bringing 
us their work [7] in our attention. This work was sup- 
ported by Grant No. 2000-2-11200-002-3 from the Basic 
Research Program of KOSEF. 



[1] Dynamics of Fractal Surfaces, edited by F. Family and 
T. Vicsek (World Scientific, Singapore, 1991); A.-L. 
Barabasi and H. E. Stanley, Fractal Concepts in Surface 
Growth (Cambridge University Press, Cambridge, 1995); 
J. Krug, Adv. Phys. 46, 139 (1997). 



[2] S. F. Edwards and D. R. Wilkinson, Proc. R. Soc. Lon- 
don, Ser. A 381, 17 (1982). 

[3] F. Family and T. Vicsek, J. Phys. A 18, L75 (1985). 

[4] Y. Kim, S. Y. Yoon, and H. Park, Phys. Rev. E 66, 
R040602 (2002). 

[5] B. D. Hughes, Random Walks and Random Environ- 
ments (Clarendon Press, Oxford, 1995); H. E. Stanley, 
K. Kang, S. Redner, and R. L. Blumberg, Phys. Rev. 
Lett. 51, 1223 (1983); 54, 1209 (1985). 

[6] M. D. Donsker and S. R. S. Varadhan, Commun. Pure 
App. Math. 32, 721 (1979). 

[7] H.-C. Jeong and J. M. Kim, Phys. Rev. E 66, 051605 
(2002). 

[8] J. M. Kim and S. Das Sarma, Phys. Rev. E 48, 2599 
(1993). 

[9] J. D. Noh, H. Park, D. Kim, and M. den Nijs, Phys. 
Rev. E 64, 046131 (2001); J. D. Noh, H. Park, and M. 
den Nijs, Phys. Rev. Lett. 84, 3891 (2000). 
[10] D.-S. Lee and M. den Nijs, Phys. Rev. E 65, 026104 
(2002). 

[11] H.-C. Jeong, J. M. Kim, H. Choi, and Y. Kim (unpub- 
lished) . 



3 



